Magnetic Polarization Currents in Double Quantum Dot Devices 
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We investigate coherent electron transport through a parallel circuit of two quantum dots, each 
of which has a single tunable energy level. Electrons tunneling via each dot from the left lead 
interfere with each other at the right lead. It is shown that due to the quantum interference of 
tunneling electrons the double quantum dot device is magnetically polarized by coherent circulation 
of electrons on the closed path through the dots and the leads. Varying the energy level of each dot 
one can make the magnetic states of the device to be either up-, non-, or down- polarization. It is 
shown that for experimentally accessible temperatures and applied biases the magnetic polarization 
currents should be sufficiently large to observe with current nanotechnology. 

PACS numbers: 73.23.Hk, 73.63. Kv, 73.40. Gk, 85.35.Ds 



The most interesting phenomena seen in mesoscopic 
electronic devices are due to the quantum coherence of 
electrons being maintained over a significant part of the 
transport process. Examples of such interference effects 
that have been observed include weak localization, 
universal conductance fluctuations, and Aharonov-Bohm 
(AB) oscillations. In 1995, Yacoby and coworkers ||] 
demonstrated the coherence of electron waves passing 
by resonant tunneling through a quantum dot (QD) in 
a double-slit type interference in a ring geometry. Re- 
cent interference experiments [^-|[ with two different 
transport paths in a ring geometry have enabled the 
realization of a phase sensitive probe of the effects of 
electron-electron interaction on the conductance oscil- 
lation such as Kondo correlations as well as the 
anomalous phase of the transmission coefficients through 
a QD [^-15|. Two QDs have also been fabricated exper- 
imentally on two different electron pathways [|l6|. These 
double quantum dot (DQD) devices provide a good op- 
portunity to test theories of resonant tunneling |]l7| , |l8| , 
cotunneling [^^|-^, and many-body correlation effects 
p2| , ^. Compared to ballistic electron interference de- 
vicespj,^, a DQD device makes it possible to manip- 
ulate the coherent tunneling of electrons through each 
dot separately by varying the gate voltages of the dots. 
Konig and Gefen |21| have discussed quantum coherence 
in DQD devices with the same energy level in each dot. 

In this Letter, we study coherent electron transport 
through two parallel QDs, each of which has a single 
tunable energy level (see Fig. |l|). Remarkably, we find 
a coherent magnetic polarization current (MFC) circu- 
lating on the closed path connecting the dots and the 
leads as a function of each dot level position. This MFC 
is induced by coherent tunneling for electron transport 
through each QD. We discuss the magnetic polarizability 
of the DQD device due to the MFC for finite temperature 
and finite applied bias. 

We start with the model Hamiltonian 



~Y£^4,d,^+Y. {ykclAa+\i.c.), (1) 



iel,2 



where Cka- and dia- are the annihilation operators with 
spin a for electrons in the leads and the dots {i — 1,2), 
respectively, ei and £2 are the level energy in each 
dot, measured, relative to the Fermi energy of the leads. 
The symmetric tunnel-coupling between the dots and 
the leads will be assumed to be independent of energy, 
1^41 = IV'I. 

The current flowing into the each quantum dot can be 
deflned as the rate of change in the number of electrons 
in a lead. At the left lead L, the total current is split 
into two local currents, Ii and l2- The commutator of 
the number operator N = X^/co-eL '-i^a-'^kcr with the Hamil- 
tonian (|^) gives rise to the current as the sum of the local 
currents through each dot. 
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with the nonequilibrium Green's function G^^ = 

i{d\^{t')cko(t))- With the Keldysh technique for non- 
linear current through the system, the local currents 
through each dot are given by P^,p7| 



= / de{fde)~fR{e))%{e), 
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where the local transmission spectral functions are de- 
fined by T,{e) = {r^G;,(£)r«G^(e)}^^ which is the i- 
th diagonal component of the matrix transmission spec- 
tral function. Here, fa{£) = /(e — /^a) is the Fermi- 
Dirac distribution function of the leads a — L, R and 
fiL — —fJ-R — eV/2 with applied bias eV between two 
leads. Due to tunneling each dot level acquires a fi- 
nite line width F = 2Tr\V\'^J\f, where TV is the density 
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of states in the leads. The matrix coupling to the leads 

is described by = = T ^ | | ^ . G;,(e) is the 

matrix dot Green's function defined in time space as 
Gl,,At-t') = -t0it-t'){{d.,{t),dlit')}). By using 
the equation of motion treatment, one can obtain the 
matrix Green's function of the dots as 



e — ei +iT 
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and G^(£) — [GJ^(£)]^. Accordingly, the local transmis- 
sion spectral functions are written by 
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Note that these can be negative. The total current is 
the sum of current through each dot I = Ii + I2 which 
is just the current conservation. This leads to the total 
transmission spectral function as T{e) = Ti{e) + T2{e), 
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We note that this is always positive. The classical ana- 
logue of our system is two resistors in parallel, h and I2 
must then both be positive. In contrast, in a quantum 
system the only constraint is that current conservation 
requires I = Ii + I2. It is not required that / > /i , l2- 
For the case of a metallic ring coupled to leads, this was 
pointed out previously by Jayannavar and Deo ||2^ 

Let us assume the cases oi I < Ii or I < I2 under 
the current conservation for /i/j < /il- For given energy 
levels (£i,£2), if /(£i,£2) < /i(£i,£2), we can assign an 
excess current lexci^ii^T)- Then we can rewrite the to- 
tal current as /(£i,£2) = /i(£i,£2) - /e2;c(£i, £2). The 
current conservation gives rise to the local excess cur- 
rent of /ea:c(ei, £2) = " ^2(ei,£2) which should circulate 
clockwise on the closed path through the dots and the 
leads. In the opposite case of / < /2, the excess current 
becomes /ea:c(£ij£2) — ~'Ii{£i,S2) circulating counter- 
clockwise on the closed path. The circulating current 
makes the device magnetically polarized. Therefore, we 
define the circulating current as a magnetic polarization 
current (MPC) Im = —lexc- We choose its direction for 
the case of / < /i as positive. It should be noted that this 
is purely a quantum coherent mesoscopic phenomena. 

Considering the transport current (TC), /, and MPC, 
Im, on an equal footing, we define the MPC as 
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with the effective spectral function, Tm(£)- The Tm{s) 
can be extracted from the following arguments. Let us re- 
call the transmission spectral functions for £1 < £2. T{e) 



has three extremum points, that is, T(£i) = T(£2) = 1 
(resonant transmission) and T{e) = (anti-resonant 
transmission), where £ — [ei -|-£2)/2. At £ = £, the anti- 
resonance of T(£) gives rise to a pronounced dip structure 
originating from the destructive interference between the 
transmissions through one QD and the other. Such an 
anti-resonant feature in a transport system with two dif- 
ferent transmission channels is well understood as the 
Fano effect Next, the two local transmission spec- 

trals of %_{£) and 72(e) have three characteristic points, 
that is, ri(£i) = 1 [r2(£2) = 1] and ri(£2) = Ti{e) = 
[72 (ei) — 72(e) = 0]. These points have nothing to do 
with resonant and anti-resonant tunneling through each 
dot. The two local transmission spectral functions only 
give us information about the local currents. Then we 
have to determine the behavior of the local spectrals in 
other energy regions. It is convenient to consider the ra- 
tio of the local transmission to the total transmission. 
The ratios are written by Ti{£)/T{e) = 1/(1 + g{£)) and 
r2(£)/r(£) = 1/(1 + g(£)-i), where gie) = [e - £i)/(£ - 
£2). For £ < £1 and e > £2, since < g{e) < 1, the ratios 
are between and 1. In these regions, there are no local 
excess currents. However, for £1 < £ < £, 7i(£)/T(£) > 1 
[r2(£)/r(£) <0] and for £ < £ < £2, Ti[e)/T{e) < 
[72(£)/T(£) > 1]. These spectral properties give rise to 
the MPC at a given energy e. Similarly, one can de- 
cide the Tm{s) for £1 > £2. Consequently, we obtain the 
7j\/(£) in terms of the local transmission spectrals as 



7'm(£) = ^(^^ - - - ^)^ (^) 

-0(£-£,)0(£--£)7-(e)}. (8) 



Equations (0) and ^ are the central result of this work. 
Note that, for a given energy level position (£i,£2) in 
each dot, the 7a/ (e) is non-zero between the two energy 
levels Ei < £ < Ej and is an anti-symmetric function with 
respect to £, Tm{£ — e) = — 7a/(£ — e)- These properties 
of Tm{e) determine the window of applied bias in which 
the MPC can be measured. 

At zero temperature, the limit of zero applied bias is 
the most simple case. The TC is proportional to the 
transmission of incoming electrons at the Fermi energy 
(£_F = 0); limey^o I = {2e/h)T{e)\e=ei.eV and the MPC 
becomes limey ^0 7m = (~2e//i)7M(£)|£=eFey. When 
the energy level of one dot is lying below the Fermi en- 
ergy and that of the other is lying above the Fermi en- 
ergy, the MPC appears to polarize the DQD device. If 
both energy levels of dots are below or above the Fermi 
energy the device is not magnetically polarized. This 
implies that the interference between the electron and 
hole channels produces the MPC. Figure ^ shows the 
magnetic polarization as a function of (£i/r,£2/r) for 
f-R < A^L- It is shown that manipulating the energy level 
position of each dot, one can magnetize the DQD devices 
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as up-, non-, or rfown-polarized. Applying a finite bias 
between the leads, the properties of Tm{s) change the 
polarization zone boundaries. A finite applied bias de- 
velops a non-polarization zone satisfying the conditions 
of -eV/2 < |ei - £2! < eV/2 or -eV/2 <ei,S2< eV/2. 
While the up-, and rfozwn-polarization zones are extended 
to the non-polarization zone of the limit of zero applied 
bias because the electron and hole channels near the 
Fermi energy within the window of the applied bias con- 
tribute to the MPC. It should be noted, when the applied 
bias is reversed to fiR > /xl, the magnetic moment of the 
device is reversed. 

To illustrate the MPC for finite temperature, we choose 
a set of energy level position (ei/F, £2/r) = (0.3,-0.9) 
which can be adjusted to other values by varying the gate 
voltages. In fact, the level positions taken in the same 
polarization zone do not affect the physics of the MPC 
but only change its amplitude. We display the MPC and 
the TC as a function of applied bias for different tem- 
peratures in Fig. ||. As the applied bias increases from 
zero bias, both the TC and the MPC increase linearly. 
The MPC is always smaller than the TC for these given 
energy level positions. However, for the case of other en- 
ergy level positions, the MPC can become larger than the 
TC (e.g., for (£i/F,e2/F) = (0.5,-0.6), Im ^51 at low 
temperatures). This linear behavior of the MPC shows 
that the MPC emerges only in nonequilibrium. Further 
increase of applied bias results in the MPC approach- 
ing its maximum value. Eventually, the disappearance 
of the MPC arises when the window of the applied bias 
becomes larger than the range of e for which Tm (e) has a 
non-zero value (-0.9 < £ < 0.3). The inset of Fig. |(b) 
shows that the TC increases non-linearly as the applied 
bias increases. In addition, compared to the MPC, the 
TC is suppressed for the bias smaller than eV = 0.6 F 
but enhanced for the bias larger than eV = 0.6 F by 
thermal effects. This originates from the fact that T{e) 
has a pronounced dip structure at £ = — 0.3F, due to 
the Fano effect. However, the anti-symmetric property 
of Tm{s) gives rise to moderate thermal suppression of 
the MPC. The thermal suppression is manifestly shown 
in the temperature dependence of the MPC and the TC 
in Fig. H The relatively large applied bias leads to the 
large amplitude of the MPC. This is consistent with the 
linear behavior of the MPC in the Im-V curve. The in- 
sets of Fig. ^(a) and (b) show, compared to the TC, the 
more rapid suppression of the MPC since Tm{s) is zero 
for low and high energies. At fc^T ~ 0.05 F, the MFCs 
at various applied biases start to be suppressed by ther- 
mal effects. For temperatures higher than ksT ~ 0.2 F, 
thermal effects wash out this novel quantum coherent 
phenomena. 

From the experimental parameters measured in Ref. 
F ~ 50fieY and A = 2.52 x lO^^m^, where A 
is a corresponding area to AB oscillation, we can esti- 
mate the amplitude of a MPC and an induced magnetic 



moment, jlu] = A ■ Im- At the point A in the inset 
of Fig. |(a), for ksT = 0.1F(T ~ 50mK), one can 
estimate Im — 0.36 nA, when eV ~ 25/^eV is applied. 
The induced magnetic moment of the device becomes 
~ 9/iB, where /^b is the Bohr magneton. Compar- 
ing this estimate of the MPC to recent measurements of 
persistent currents |25j suggests that the effects we are 
discussing can be observed with existing nanotechnology. 

In summary, we studied coherent electron transport 
through two parallel quantum dots, each of which has 
a single tunable energy level. By changing these energy 
levels in the DQD device one can vary the sign and mag- 
nitude of the magnetic polarization current induced by 
quantum interference effects. This current is sufficiently 
large that it should be experimentally observable. 
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[1] See, for a review, Y. Imry, Introduction to Mesoscopic 
Physics (Oxford University Press, Oxford, New York 
1997); S. Datta, Electronic Transport in Mesoscopic Sys- 
tems (Cambridge University Press, New York 1995). 

[2] A. Yacoby et ai, Phys. Rev. Lett. 74, 4047 (1995). 

[3] R. Schuster et al, Nature 385, 417 (1997). 

[4] W.G. van der Wiel et al, Science 289, 2105 (2000). 

[5] Y. Ji et al. Science 290, 779 (2000). 

[6] U. Gerland et al, Phys. Rev. Lett. 84, 3710 (2000). 

[7] W. Hofstetter, J. Konig, and H. Schoeller, Phys. Rev. 
Lett. 87, 156803 (2001). 

[8] K. Rang, and S.-C. Shin, Phys. Rev. Lett. 85, 5619 
(2000). 

[9] S. Y. Cho et al, Phys. Rev. B 64, 033314 (2001). 
[10] A. L. Yeyati, and M. Biittiker, Phys. Rev. B 52, R14360 
(1995). 

[11] G. Hackenbroich, and H.A. Weidenmiiller, Phys. Rev. 

Lett. 76, 110 (1996). 
[12] C. Bruder, R. Fazio, and H. Schoeller, Phys. Rev. Lett. 

76, 114 (1996). 

[13] C.-M. Ryu and S. Y. Cho, Phys. Rev. B 58, 3572 (1998). 
[14] H.-W. Lee, Phys. Rev. Lett. 82, 2358 (1999). 
[15] R. Baltin and Y. Gefen, Phys. Rev. Lett. 83, 5094 (1999). 
[16] A.W.HoUeitner et al, Phys. Rev. Lett. 87, 256802(2001). 
[17] T. V. Shahbazyan, and M.E. Raikh, Phys. Rev. B 49, 
17123 (1994). 

[18 
[19 
[20 

[21 
[22; 



B. Kubala and J. Konig, :ond-mat/0110524 



H. Akera, Phys. Rev. B 47, 6835 (1993). 
D. Loss, and E.V. Sukhorukov, Phys. Rev. Lett. 84, 1035 
(2000). 

J. Konig, and Y. Gefen, Phys. Rev. Lett. 86, 3855 (2001). 
D. Boese, W. Hofstetter, and H. Schoeller, Phys. Rev. B 
64, 125309 (2001). 



[23] W. Hofstetter and H. Schoeller, ;ond-mat/0108359 



3 



[24] D. Mailly, C. Chapelier and A. Benoit, Phys. Rev. Lett. 

70 2020 (1993). 
[25] W. Rabaud et ai, Phys. Rev. Lett. 86, 3124 (2001). 
[26] Y. Meir, and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 

(1992). 

[27] S. Y. Cho, K. Rang and C.-M. Ryu, Phys. Rev. B 60, 
16874 (1999). 

[28] A. M. Jayannavar, and P. Singha Deo, Phys. Rev. B 51, 

10175 (1995). 
[29] U. Fano, Phys. Rev. 124, 1866 (1961). 

QD 1 





£2i 




QD2 



FIG. 1. A double quantum dot (DQD) device. Both dots 
are tunnel-coupled to the left and right leads. The leads are 
characterized by the chemical potentials, jj,L and fiR. The 
tunneling amplitudes between the dots and the leads are de- 
noted by r. The energy level position in each dot is measured 
as ei and £2 from the Fermi energy in the leads. 




FIG. 2. Magnetic polarization of the double quantum dot 
device as a function of the energy level positions of each dot 
(£i/r,£2/r), in the limit of zero applied bias at zero tem- 
perature for < fiL- The vertical arrows stand for the 
magnetic moment of the DQD device whose length and direc- 
tion depends on the amplitude and direction of the magnetic 
polarization current, respectively. 




FIG. 3. (a) Magnetic polarization current and (b) its ra- 
tio to the transport current as a function of applied bias at 
(£i/r,e2/r) = (0.3,-0.9) for various temperatures ksT. In 
the insets, the currents are shown in a linear scale. 
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FIG. 4. Temperature dependence of (a) the magnetic po- 
larization current and (b) its ratio to the transport current at 

(ei/r,£2/r) — (0.3,-0.9) for different values of the applied 
bias eV. In the insets, the currents are shown in a linear scale. 
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